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The relativistic field theory model of the deuteron (RFMD) is applied to the cal- 
qq \ culation of the astrophysical factor 5 , pp (0) for the process of the solar proton burning 

p + p— >D + e + + f e and the cross sections for the disintegration of the deuteron by 
photons 7 + D — > n + p and anti-neutrinos v e + D — ► e + + n + n. Our theoretical 
value of the astrophysical factor S , pp (0) = 4.02 x 10 -25 MeVb agrees with the clas- 
sical result obtained by Bahcall and Kamionkowski 5^(0) = 3.89 x 10~ 25 MeVb in 
the potential model approach (PMA). The cross sections for the disintegration of 
the deuteron by photons and anti-neutrinos calculated near thresholds are in good 
agreement with the PMA. An extrapolation of the cross sections for energies far 
^ . from thresholds is suggested and related to the inclusion of form factors describing 

spatial smearing of the deuteron and the NN system. The extrapolated cross section 
for the disintegration of the deuteron by anti-neutrinos agrees with that calculated 
in the PMA in the anti-neutrino energy region from threshold up to Ep c = lOMeV. 
The extrapolated cross section averaged over the reactor anti-neutrino energy spec- 
trum is obtained in agreement with the experimental data. It is shown that the 
RFMD enables to describe elastic low-energy NN scattering in accordance with 
low-energy nuclear phenomenology. 
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1 Introduction 



The relativistic field theory model of the deuteron (RFMD) formulated in Refs. [1,2] 
has been applied to the calculation of the reaction rate of the neutron-proton radiative 
capture n + p — > D + 7 and the astrophysical factor 5 , pp (0) of the solar proton burning 
p + p — > D + e + + v c [2]. Some mistakes, which had been made for the first calculation 
[2], have been then partly, mainly for the reaction rate of the neutron-proton radiative 
capture, corrected in Ref. [3]. However, our conclusion concerning the value of the astro- 
physical factor S'pp(O) [3] is still erroneous. 

In this paper we would like to amend our results obtained in Refs. [1-3] and to apply 
the RFMD to the calculation of the cross sections for the disintegration of the deuteron 
by photons 7 + D — > n + p and anti-neutrinos v c + D — > e + + n + n. The reaction 
of the disintegration of the deuteron by anti-neutrinos u c + D — > c + + n + n is caused 
by the charged weak current and valued, in the sense of charge independence of the 
weakinteraction strength, to be equivalent to the observation of the reaction of the solar 
proton burning p + p — > D + e + + v c in the terrestrial laboratories [4]. Experimentally 
the reaction of the disintegration of the deuteron by anti-neutrinos v e + D — > e + + n 
+ n induces itself by reactor anti-neutrinos with an equilibrium energy spectrum [5,6]. 
Therefore, experimental data on the reaction z/ c + D — >e + + n + n are given in the form 
of the cross section averaged over the reactor anti-neutrino energy spectrum [4,7-9]. 

We show that all processes under consideration can be described in the RFMD in 
agreement with the potential model approach (PMA) in spite of completely different dy- 
namics of strong low-energy nuclear interactions. Indeed, in the RFMD [1,2] the physical 
deuteron appears through long- wavelength vacuum fluctuations of the proton and the 
neutron field in the one-nucleon loop approximation. In terms of one-nucleon loop ex- 
changes we describe in the RFMD a non-trivial wave function of the relative movement 
of the nucleons inside the physical deuteron. Therefore, the physical deuteron couples to 
nucleons and other particles only through one-nucleon loop exchanges. 

In order to couple to the deuteron through the one-nucleon loop exchange the nucleons 
should pass through intermediate interactions providing low-energy transitions N + N — > 
N + N. In a quantum field theory approach such interactions should be induced by meson 
exchanges. As the nucleons couple at low energies, the main contribution should come 
from the one-pion exchange. The contributions of heavier meson exchanges can be taken 
into account effectively by integrating them out. 

Since in the reactions p + p^D + e + + z/ c , n + p^D + 7, 7+ D — > n + p and u c 
+ D^e + + n + n the nucleons couple in the 1 So-state, the low-energy transitions N + 
N — > N + N can be described by the effective local four-nucleon interactions [1,2]: 

£ NN^NN (x) = G ^ mx)l ^{x)mx)ri b n{ X )] 

+- [n(x)^^ b n c (x)][n c (x)Y^n(x)] + - [p(x)^^ b p c (x)][p c (x)^^p(x)} 

+ (7 M 7 5 ®7 M 7 5 ^7 5 ®7 5 )}, (1-1) 

where n(x) and p(x) are the operators of the neutron and the proton interpolating fields, 
n c (x) = Cn T (x), etc., then C is a charge conjugation matrix and T is a transposition. 
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The effective coupling constant Gvnn is defined by 

G.nn = fjjf - ^ = 3.27 x 10- 3 MeV" 2 , (1.2) 
4il^ M N 

where (JVnn = 13.4 is the coupling constant of the 7rNN interaction, M n = 135 MeV is 
the pion mass, M p = M n = Mn = 940 MeV is the mass of the proton and the neutron 
neglecting the electromagnetic mass difference, which is taken into account only for the 
calculation of the phase volumes of the final states of the reactions p + p^D + e + + 
v e and u e + D — > e + + n + n, and a np = (—23.748 ± 0.010) fm is the S-wave scattering 
length of the np scattering in the 1 So-state [10]. 

The first term in the effective coupling constant G^nn comes from the one-pion ex- 
change for the squared momenta transfer — q 2 much less than the squared pion mass 
— q 2 <C M 2 and the subsequent Fierz transformation of the nucleon fields [1,2]. We 
should emphasize that due to Fierz transformation the effective NN interaction caused 
by the one-pion exchange contains a few contributions with different spinorial structure, 
we have takeninto account only those terms which contribute to the 1 So-state of the NN 
system. The second term is a phenomenological one representing a collective contribution 
caused by the integration over heavier meson fields like scalar mesons <r(700), ao(980) and 
/o(980), vector mesons p(770) and a; (780) and so on. This term is taken in the form used 
in the Effective Field Theory (EFT) approach [11-13]. The effective interaction Eq. (|1 . 1| ) 
is written in the isotopically invariant form, and the coupling constant G^nn can be never 
equal zero at a np ^ due to negative value of a np imposed by nuclear forces, i.e., a np < 
[14]. 

In the low-energy limit the effective local four-nucleon interaction Eq. fll.lj) vanishes 
due to the reduction 

[iV(x)7 M 7 5 iV c (x)][iV c (x)7V^(x)] -> -[N(x)-f 5 N c (x)][N c (x)-f 5 N(x)], (1.3) 

where N(x) is the neutron or the proton interpolating field. Such a vanishing of the 
one-pion exchange contribution to the NN potential is well-known in the EFT approach 
[11-13] and the PMA [14]. In power counting [11-13] the interaction induced by the one- 
pion exchange is of order 0(k 2 ), where k is a relative momentum of the NN system. The 
former is due the Dirac matrix 7 s which leads to the interaction between small components 
of the Dirac bispinors of the nucleon wave functions. 

Thus, if either in the PMA or the EFT approach the effective local four-nucleon 
interaction Eq. ( |1.1|) would be applied to the description of the deuteron coupled to the 
nucleons, the contribution would be scarcely significant. Therefore, both in the PMA 
and the EFT for the correct description of strong low-energy nuclear forces one needs to 
include an effective phenomenological NN potential, for instance, the Argonne t>ig [15]. 
The one-pion exchange contribution is considered as a perturbation. 

In the RFMD due to the one-nucleon loop exchange the contributions of the inter- 
actions [A(x)7 Ai 7 5 A^ c (a;)][A c (x)7 M 7 5 iV(x)] and [A(x)7 5 A^ c (x)][A c (x)7 5 A^(2;)], or shortly 
7^7 5 ® 7^7 5 and 7 5 ® 7 5 , to the amplitudes of nuclear processes are different and do 
not cancel each other in the low-energy limit. This is completely a peculiarity of one- 
nucleon loop diagrams related to one-fermion loop anomalies [16,17,2]. For instance, in 
the case of the neutron-proton radiative capture and the photomagnetic disintegration of 
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the deuteron the amplitudes of the processes are defined by the triangle one-nucleon loop 
diagrams with AVV (axial-vector-vector) and PVV (pseudoscalar-vector-vector) vertices 
[2] caused by 7 M 7 5 <E> 7 M 7 5 and 7 5 <g> 7 5 interactions, respectively. These diagrams are well- 
known in particle physics in connection with the Adler-Bell-Jackiw axial anomaly [16] 
which plays a dominant role for the processes of the decays 7r° — > 77, uo — > 7r°7 and so on 
[16]. The results of the calculation of these diagrams differ each other. Hence, they give 
different contributions to the amplitudes of the processes and do not cancel themselves in 
the low-energy limit. 

Then, the amplitudes of the solar proton burning and the anti-neutrino disintegration 
of the deuteron are defined by the one-nucleon loop diagrams with AAV and APV vertices 
caused by 7^7 5 <8> 7 M 7 5 and 7 s <g> 7 s interactions, respectively. The contribution of the 
diagrams with APV vertices turns out to be divergent and, therefore, negligibly small 
compared with the contribution of the diagrams with AAV vertices [2], which contains 
non-trivial convergent part related to one-fermion loop anomalies [2,17]. 

As a result in the low-energy limit amplitudes of nuclear processes described by the 
RFMD contain only large components of Dirac bispinors of wave functions of nucleons. 
This provides an effective enhancement of the one-pion exchange for the description of 
strong low-energy interactions of the nucleons in the 1 S -state. Thus, in the RFMD the 
dynamics of strong low-energy nuclear interactions caused by one-nucleon loop exchanges 
is the point of a dominant role of the one-pion exchange contribution to the effective low- 
energy interactions of nucleons in the 1 So-state. 

The paper is organized as follows. In Sect. 2 we calculate the astrophysical factor for 
the solar proton burning. In Sect. 3 we calculate the cross section for the photomagnetic 
disintegration of the deuteron near threshold and formulate the extrapolation procedure 
for the cross section for energies far from threshold. In Sect. 4 we compute the cross section 
for the anti-neutrino disintegration of the deuteron near threshold, extrapolate this cross 
section for energies far from threshold and average the extrapolated cross section over the 
reactor anti-neutrino energy spectrum. In Conclusion we discuss the obtained results and 
a justification of the RFMD. In Appendix A we adduce the detailed calculation of the 
amplitude of the solar proton burning. In Appendix B we show that in the RFMD one 
can describe low-energy elastic NN scattering with non-zero effective range in accordance 
with low-energy nuclear phenomenology. 



2 Astrophysical factor for the solar proton burning 

In the RFMD the amplitude of the solar proton burning is defined by one-nucleon loop 
diagrams [2]. The detailed calculation of the amplitude of the solar proton burning is 
given in Appendix A. The result reads 

iM(p + p -> D + e+ + v e ) = C(rj) G Y g A M N G wNN ^ 

x e;(k D ) [u{K e )^{\ - 7 5 MM] [u^tVpi)], (2.1) 

where Gy = Gf cos $c with G F = 1.166 x 1(T 5 GeV~ 2 and tf c are is the Fermi weak 
coupling constant and the Cabibbo angle cos-# c = 0.975. Then g A = 1.260 ± 0.012 
describes the renormalization of the weak axial hadron current by strong interactions 
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[10], gv is the phenomenological coupling constant of the RFMD related to the electric 
quadrupole moment of the deuteron: g Y = 2ir 2 Q I )M§ l [2] with Q D = 0.286 fm 2 [10], e*(/c D ) 
is a 4-vector of a polarization of the deuteron and u(k Uc ), v(k e +), u c (p 2 ) and u{pi) are 
the Dirac bispinors of neutrino, positron, and two protons, respectively. For the binding 
energy of the deuteron we use the value Ed = 2.225 MeV [10]. The Coulomb repulsion 
between protons is taken into account only in terms of the Gamow penetration factor [2,3] 
C{rf) = a/27t?7 exp(— irrf) depending on the relative velocity of the protons v as 77 = a/v 
and a = 1/137 is the fine structure constant. 

The cross section for the low-energy p + p— > D + e + + v e reaction is defined 

a(pp - De+z/ e ) = - — I— f |^( p + p ^D + e+ + z/ e )| 2 

x(2^S^ + K, + K.- Pl - P2) ^_^ E -^-, (2.2) 

where v is a relative velocity of the protons and Ei(i — 1,2) are the energies of the protons 
in the center of mass frame. 

Then, \A4(p + p — ► D + e + + i/ e )| 2 is the squared amplitude averaged over polariza- 
tions of protons and summed over polarizations of final particles: 

|.M(p + P-D + e+ + ^ = C\rj) G Y g 2 A M* G; NN ^ x ( - f + ^ 

xtr{(-m e + k e +ha(l ~ 7 5 )^ c 7/3(l - 7 5 )} x \ x tr{(M N - p 2 ) 7 5 (M N + p^ 5 }, (2.3) 

where m e = 0.511 MeV is the mass of positron, and we have used the relation gy/n 2 = 
2Q D M 2 [2}. 

In the low-energy limit the computation of the traces yields 



- 9^ + ^ ] x tr{(-m e + k e+ )l a {l ~ l 5 )k»M± ~ 7 5 )} 



= 24 \E e+ E„ e -^k e+ -k„ A 

~ x tr{(M N -p 2 ) 7 5 (M N +p 1 ) 7 5 } = 2M 2 , (2.4) 

where we have neglected the relative kinetic energy of the protons with respect to the 
mass of the proton. 

Substituting Eq. in Eq. (pj) we get 



\M{ V + P - D + e+ + u e )\ 2 = C 2 ( V ) Gy g 2 A M* G; NN — ^ ( E e+ E„ e - -k e+ -k Vc ). (2.5) 



The integration over the phase volume of the final De + i/ e -state we perform in the non- 
relativistic limit 

/ W pfe (ife ^ S * Uk ° + ~ » " ^ ~ ^ ■ l - ) 

1 rW+T pp , fW 1 ^5 



32tT 3 M, 



where W = e B - (M n - M p ) = (2.225 - 1.293) MeV = 0.932 MeV, T pp = M N t» 2 /4 is 
the kinetic energy of the relative movement of the protons, and £ = m e / (W + T pp ). The 
function /(£) is defined by the integral 



/(0 = 30 jf 1 y^ 2 - £ 2 a: (1 - xfdx = (1 - - 4£ 4 ) 0^ 

= 0.222 (2.7) 



15 , 4 „ /1 + V / W 



and normalized to unity at £ = 0. 

Thus, the cross section for the solar proton burning is given by 

a(pp -> De + z/ e ) = — a — — G wNN + T pp ) / 



v 2 320 vr 4 ^ v pp/ J V W + T t 

^pp(^pp) — 27T77 



PP/ 



e-^ /. (2.8) 
J pp 

The astrophysical factor S PP (T PP ) reads 

Spp(Tpp) = a Mgg^ G J NN {w + Tpp)5 , , „ 9) 

At zero kinetic energy of the protons T pp = the astrophysical factor S , pp (0) is given by 
9g\ G 2 Q D M± n2 5 Jm e \_ Ano ^. n _ 25 



5pp(0) = a i28b^ G ' NN w f [w) =402 x 10 MeVb - ( 2 - 10 ) 

The value S pp (0) = 4.02 x 10~ 25 MeVb agrees with the value S pp (0) = 3.89 x 10~ 25 MeVb 
obtained by Kamionkowski and Bahcall in the PMA [18] and S pp (0) = 4.05 x 10" 25 MeVb 
having been calculated recently in the EFT approach [19]. 

Since due to charge independence of the weak interaction strength the reaction of 
the anti-neutrino disintegration of the deuteron u e + D — ► e + + n + n is valued as a 
terrestrial equivalent of the solar proton burning, for the justification of the validity of 
our result Eq. ( |2.10| ) we suggest to calculate the cross section for the disintegration of the 
deuteron by reactor anti-neutrinos u c + D — > e + + n + n. This calculation is carried 
out in Sect. 4. However, the RFMD with the effective local four-nucleon interaction 
Eq. (11. If) allows to calculate the cross section for the reaction z/ e + D — >e + + n + n 
for anti-neutrino energies close to threshold. As it is shown in Sect. 4 this cross section 
agrees well with the result obtained in the PMA. For the comparison of our result with 
the experimental data represented in the form of the cross section averaged over the 
anti-neutrino energy spectrum we need to extrapolate the cross section calculated near 
threshold for the anti-neutrino energies far from threshold. In order to formulate such an 
extrapolation procedure we turn to the consideration of the photomagnetic disintegration 
of the deuteron 7 + D — > n + p. 
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3 Photomagnetic disintegration of the deuteron 



First, let us consider the cross section for the photomagnetic disintegration of the deuteron 
7 + D — > n + p near threshold. The cross section can be easily obtained by using the 
results of Ref. [2,3]. It reads 

ct^(uj) = a (|-J kr D , (3.1) 

where k = ^Jm^(uj — e-£>) is the relative momentum of the np system, uj is the energy of 
the photon and td = lfy/e^M^ = 4.315 fm is the radius of the deuteron, and a is given 
by ' 

a ° = { ^ - w g ' nn4/2m ^ 2 = 5 - 9 x 10 " 7cm2 - (3 ' 2) 

The cross section a lD (u), calculated in the PMA near threshold has the same form as 
Eq. ( |3.1|) but with uq given by [14] 



0o 



2ira 
3M§ 



(Jh - ^) 2 (l - Ahp^dMn) 2 = 6.3 x 10- 27 cm 2 . (3.3) 



It is seen that o"o defined by Eq. (|3.2|) and Eq. ( |3.3|) agree within an accuracy better than 
10%. The analogous agreement can be drawn out from the comparison of the reaction 
rates for the neutron-proton radiative capture for thermal neutrons [3,14]: 

2 25aQ D 2 3 ln -30„_2 



(/i p - ^Y—r Y- Gi NN 4m n = 2.0 x 10 



cm 



(/^-^Vl 1 "^^) (it^J = 2 -2 x 10~ 30 cm 2 , 



where v and k are a relative velocity and a relative 3-momentum of the np system. Thus, 
near threshold the RFMD supplemented by the effective local four-nucleon interaction 
Eq. provides a dynamics of strong low-energy nuclear forces describing the cross 

sections for the neutron-proton radiative capture and the photomagnetic disintegration 
of the deuteron in good agreement with the PMA. 

Now let us proceed to the formulation of the extrapolation procedure. For this aim we 
suggest to consider the cross section for the photomagnetic disintegration of the deuteron 
calculated in the PMA for the photon energies far from threshold [14] 

= it^i^ = ff °fe) kr " F ^ (3 ' 5) 

where the function F® p (k 2 ) is defined as 

FR(k 2 ) = Kt^ (3-6) 

npV ; 1 + r^k 2 1 + a 2 Qp k 2 v ; 

The cross section Eq. (|3.5| ) differs from the cross section Eq. ( |3.1|) calculated near threshold 
by the factor F^ p (k 2 ), which can be considered as a form factor taking into account a 
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spatial smearing of the deuteron through the factor 1/(1 + r^k 2 ) and the np system 
through the factor 1/(1 + al p k 2 ). 

In order not to come into the contradiction with the PMA we suggest to extrapolate 
the cross section Eq. ( |3.1| ) for the photon energies far from threshold by means of the form 
factor F^ p {k 2 ). This extrapolation assumes that the cross section for the photomagnetic 
disintegration of the deuteron calculated in the RFMD has the form Eq. ( |3.5|) with cto 
defined by Eq. (|3.2| ). Such an extrapolation applied to the cross section for the anti- 
neutrino disintegration of the deuteron should lead to the appearance of the form factor 

F ^ = TTTgFTW (37) 

describing a spatial smearing of the deuteron through the factor 1/(1 + r^k 2 ) and the nn 
system through the factor 1/(1 + a 2 n k 2 ), where a nn = — 17fm is the S-wave scattering 
length of the nn scattering in the 1 S -state [10]. 

We should accentuate that we have discussed the photomagnetic disintegration of the 
deuteron in order to draw out the hint for the formulation of the extrapolation procedure 
for the cross section for the disintegration of the deuteron by anti-neutrinos. We are not 
aiming here to compute a total cross section for the disintegration of the deuteron by 
photons. In fact, it is well known [14] that a photomagnetic part predominates only at 
small relative momenta of the np pair at krv <C 1 and the photoelectric part becomes 
important at krj) > 1, i.e., uj > = 4.45 MeV. 

Also the computation of the reaction rate for the neutron-proton radiative capture 



given by Eq. Q3.4Q can be considered as a lowest approximation. Certainly, the cross 
section a np (k) = 276 mb computed for thermal neutrons at laboratory velocities v/c = 
7.34 • 10 _ 6 (the absolute value is v = 2.2 • 10 5 cm/sec) [3] agrees reasonably well with 
the experimental data cr^ p {k) = (334.2 ± 0.5) mb [20]. However, the central theoretical 
value is 20% smaller than the central experimental value. Such a problem has been 
solved for the first time within the PMA by Riska and Brown [21] who showed that 
the discrepancy of order 10% between the theoretical value of the cross section a np (k) = 
(302.5±4.0) mb [14] calculated in the PMA and the experimental value a^ p (k) = (334.2± 
0.5) mb can be explained by exchange-current contributions. In the EFT approach and 
Chiral perturbation theory the same result has been obtained by Park et al. [22]. The 
investigation of such fine effects of the neutron-proton radiative capture and incorporation 
of Chiral perturbation theory into the RFMD is in the programme of further applications 
of the RFMD to the processes of low-energy interactions of the deuteron. 



4 Anti— neutrino disintegration of the deuteron 

The calculation of the amplitude of the anti-neutrino disintegration of the deuteron is 
analogous to the calculation of the amplitude of the solar proton burning. The effective 
Lagrangian responsible for the transition z/ e + D — > e + + n + n is given by Eq. (A. 49) 

^p cD ^o+nn(a;) = gAG n ^—^=—^D^ u (x) [n(x)^ 5 n c (x)) [ip Uc (x)j u (1 - -f 5 )tp e (x)}. 
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The amplitude of the v e + D — > e + + n + n process reads 

iM{v c + D -> e + + n + n) = - g A M N G wNN —j= — 

xe M (Q) [tJ(A*,)r(l - 7>(M] [u( Pl ) 7 V(p 2 )], (4.1) 

where v(k Pe ), v(k e +), u{p\) and u c (p2) are the Dirac bispinors of the anti-neutrino, positron 
and neutrons, e^(Q) is the 4-vector of the polarization of the deuteron. We have taken 
into account that u(p 2 ) , y 5 u c (pi) = —u(pi) r y 5 u c (p2)- The amplitude Eq. ( |4.1| ) squared, 
averaged over polarizations of the deuteron and summed over polarizations of the final 
particles reads 

l-M(^ + D^e+ + n + n)P = g 2 A M^Gl NN 144( ^ D ( Ee+ E 9e - l -k e+ • %\ . (4.2) 

Due to charge independence of the weak interaction strength the matrix element Eq. 
is related to the matrix element of the solar proton burning Eq. ( j2.5|) by the relation 

1 

3 " C\ri) 



\M(is c + D ^e+ + n + n)| 2 = - x x \M(p + p -> D + e+ + z/ e )| 2 , (4.3) 



where 8/3 is a combinatorial factor. This relation means that the dynamics of strong 
low-energy nuclear forces governing the processes of the anti-neutrino disintegration of 
the deuteron z/ e + D^n + n + e + and the solar proton burning p + p— >D + e + + z/ c 
should have the same origin. 

The expression Eq. (|4.3|) extrapolated for energies far from threshold according to the 
procedure suggested in Sect. 3 reads 

= 9lM 6 A^^[E c+ E Dc - \k e+ ■ (4.4) 



\M(p e + D -»• e+ + n + n) 



The form factor F^ n (k 2 ), describing a spatial smearing of the deuteron and the nn system, 
is given by Eq. ( |3.7|) 

C(* 2 : 1 1 



l + r 2 P l + a 2 nn k 2 ' 

where k = \/ M^T nn is the relative momentum and T nn is the kinetic energy of the relative 
movement of the nn system, and a nn = — 17 fm [10] is the S-wave scattering length of the 
low-energy elastic nn scattering in the 1 So-state. A much more complicated extrapolation 
of the form factors of the amplitude of the disintegration of the deuteron by anti-neutrinos 
has been suggested by Mintz [23]. 

The cross section for the process u e + D — > e + + n + n is defined by 



cr DeD (E* e ) = 4^-; J |A4(P c + D^e+ + n + n)P 

2 ( ) * ( \Q + K - Pl -P2- K + ) (27r)32Ei (27r)3 ^ 2 {27T)32Ee+ , (4.5) 
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where Ed, E Ve , Ei, E 2 and E e + are the energies of the deuteron, the anti-neutrino, the 
neutrons and the positron. The integration over the phase volume of the (nne + )-state we 
perform in the non-relativistic limit and in the rest frame of the deuteron, 



1 f d 3 Pl d 3 p 2 d 3 k c+ , 4 (4) , 

2 J (2^)32^ (2vr)32E 2 (2 7 r)32E e+ (27r) 6 {Q + ^ ~ Pl ~ P2 ~ h 




T P + — T„ 



1024vr 



(4.6) 



where T e + and m e = 0.511 MeV are the kinetic energy and the mass of the positron, E t ^ 
is the anti-neutrino energy threshold of the reaction z/ e + D— >e + + n + n and is given 
by E th = e D + m e + (M n -M p ) = (2.225 + 0.511 + 1.293) MeV = 4.029 MeV. The function 
f(y), where y = E V J E th , is defined as 

8 ) , 

f(y) = - dxy/x(l-x) FZ(M N E th (y - 1) x) 

IT J v 





x f 1 + — (y - 1)(1 - x)) Jl + ^(y - 1)(1 - x), (4.7) 



where we have changed the variable T nn = (E Pc — E t ^) x. The function f(y) is normalized 
to unity at y — 1, i.e., at threshold E Pc = E t ^. Thus, the cross section for the anti-neutrino 
disintegration of the deuteron reads 

a^(E Pc ) = a (y-l) 2 f(y), (4.8) 

where a is defined by 
-0 = Qd Gl NN ( ( = (4-53 ± 0.86) x lO^ 3 cm 2 . (4.9) 



512vr 4 \M N J \E, 



th 



Here ±0.86 describes the assumed theoretical uncertainty of the RFMD which is about 
19 %[]. The value a = (4.53 ±0.86) x 10" 43 cm 2 agrees with the value a = (4.68 ±1.14) x 
10" 43 cm 2 obtained in the PMA [24,25] (see Fig. 7 of Ref. [4]). 

The experimental data on the anti-neutrino disintegration of the deuteron are given 
in terms of the cross section averaged over the reactor anti-neutrino energy spectrum 
per anti-neutrino fission in the energy region of anti-neutrinos E^ < E Pc < 10 MeV: 
< a^ D (E Dc ) > cxp = (1.5 ±0.4) x 10- 45 cm 2 /z/ e fission [7], < a^ D (E Pc ) > cxp = "(0.9 ± 0.4) x 
10~ 45 cm 2 /^ fission [8] and < o^ D (E Vc ) > cxp = (1.84 ± 0.04) x 10~ 45 cm 2 /z/ e fission [9]. 

1 This is the improved estimate of the theoretical uncertainty. The former was 30% [3], 
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The cross section < a UcD (E Pc ) >, calculated in the RFMD, extrapolated and averaged 
over the reactor anti-neutrino Avignone-Greenwood energy spectrum [5,6] in the energy 
region E th < E Pc < 10 MeV, is given by 

2.482 

<a^ D (E Pc )> = "7t f dye- b y a (y-l) 2 f(y) = 

= (2.10 ± 0.40) x 1(T 45 cm 2 / P e fission, (4.10) 

where a = 17.8 E th = 71.72, b = 1.01 E th = 4.07, and N Pc = 6 is the number of anti- 
neutrinos per fission [5,6]. The theoretical value Eq. ( [4 . 1 Q| ) agrees well with the experi- 
mental values given by Reines et al. [7] and Russian groups [9], while the agreement with 
the value given by Reines et al. [8] is only qualitative. 

Thus, the calculation of the cross section for the anti-neutrino disintegration of the 
deuteron in agreement with the PMA and the experimental data confirms our result 
obtained for the astrophysical factor 5 , pp (0) = 4.02 x 10 _25 MeVb for the solar proton 
burning. 



5 Conclusion 

We have shown that the RFMD supplemented by the local four-nucleon interaction 
Eq. ( [1.1|) describes the processes of the solar proton burning p + p — > D + e + + u c , 
the neutron-proton radiative capture n + p — > D + 7, the disintegration of the deuteron 
by photons 7 + D — > n + p and anti-neutrinos u c + D — > e + + n + n in good agree- 
ment with the PMA. The astrophysical factor S pp (0) = 4.02 x 10 _25 MeVb for the solar 
proton burning, the reaction rates for the neutron-proton radiative capture and the cross 
sections for the disintegration of the deuteron by photons and anti-neutrinos calculated 
near thresholds of the reactions agree with the results obtained in the PMA within an 
accuracy better than 10%. 

In order to compare our result for the cross section for the disintegration of deuteron 
by anti-neutrinos with experimental data we suggested the procedure of the extrapolation 
of the cross section calculated near threshold to the energy region far from threshold. Of 
course, such an extrapolation is not unique. Therefore, we have suggested to formulate 
the extrapolation procedure fitting the cross section for the photomagnetic disintegration 
of the deuteron calculated in the PMA. This extrapolation assumes the multiplication of 
the cross section, calculated in the RFMD near threshold, by the form factor 

F ™ (fc2) = i+^pi+4 n p (5 - i} 

describing a spatial smearing of the deuteron by a factor 1/(1 +r 2 ) /c 2 ) and the NN system 
by a factor 1/(1 + a^ N fc 2 ), where k is a relative momentum of the NN system. A much 
more complicated extrapolation of the amplitude of the process z/ e + D^e + + n-|-n 
has been suggested by Mintz [23]. 

The extrapolated cross section for the disintegration of the deuteron by antineutrinos 
averaged over the reactor anti-neutrino energy spectrum < a UeD (E Pe ) >= (2.10 ±0.40) x 
10 _45 cm 2 / v e fission agrees well with the experimental data [7,9]. 
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The cross section for the disintegration of the deuteron by anti-neutrinos v c + D — > 
c + + n + n in dependence on the anti-neutrino energy has been calculated recently in 
the PMA in Ref. [26]. Since near threshold our cross section agrees well with the PMA 
result, for the verification of the extrapolation procedure we can, say, compare the cross 
section at E Pc = 10 MeV: 

v PcD (E Pc )\ Epc = 1QMeV = 1.02 x 10~ 42 cm 2 , RFMD, 

* PcD (E Pc )\ E 1QMeV = 1.13 xHT 42 cm 2 , PMA. (5.2) 

The visible agreement with the PMA result can serve too as a confirmation of a validity 
of our extrapolation procedure. 

The RFMD describing a relative movement of the nucleons inside the deuteron in terms 
of one-nucleon loop exchanges suggests a dynamics of strong low-energy nuclear forces 
completely different to the PMA and the EFT approach. The neutron-proton-deuteron 
vertices are point-like and defined by a phenomenological local conserving nucleon current 
J^(x) = —igy[p(x) , ~f fJ "n c (x) —n(x)'~f fl p c (x)} , i.e., d^J^{x) = 0, accounting for spinorial and 
isotopical properties of the deuteron, and gy is a dimensionless phenomenological coupling 
constant. The deuteron is represented by a local field operator D^(x) (or D' (x)), the 
action of which on a vacuum state annihilates (or creates) the deuteron. The low-energy 
parameters of the deuteron such as the binding energy e D , the electric quadrupole Q D and 
the anomalous magnetic dipole /«d moments are induced by vacuum fluctuations of the 
neutron and the proton fields in a quantum field theory way. The description of strong 
low-energy nuclear forces in terms of one-nucleon loop exchanges provides opportunity 
to convey in nuclear physics a huge experience of fermion anomalies [16,17,27] which had 
been stored in particle physics from the paper by Adler [16] concerning the derivation 
of the anomalous contribution to the axial Ward identity. In the area of low-energy 
interactions of low-lying mesons such an experience has been focused mainly upon the 
derivation of anomalous contributions [28-30] to Effective Chiral Lagrangians [31]. 

The main problem which we encounter for the practical realization of the derivation 
of effective Lagrangians of low-energy interactions of the deuteron through one-nucleon 
loop exchanges lies in the necessity to satisfy requirement of locality of these interactions 
related to the condition of microscopic causality in a quantum field theory approach [32]. 
Since in the RFMD one-nucleon loop diagrams are defined by the point-like vertices 
and the Green functions of free virtual nucleons with constant masses, there is only a 
naive way to satisfy requirement of locality of effective interactions through the formal 
application of the long-wavelength approximation to the computation of one-nucleon 
loop diagrams. This approximation implies the expansion of one-nucleon loop diagrams 
in powers of external momenta by keeping only the leading terms of the expansion. Of 
course, the application of such anapproximation to the computation of one-nucleon loop 
diagrams, when on-mass shell the energy of the deuteron exceeds twice the masses of 
virtual nucleons, can seem rather unjustified. 

However, in this connection we would like to recall that the analogous problem en- 
counters itself for the derivation of Effective Chiral Lagrangians [28,31] within effective 
quark models motivated by QCD like the extended Nambu-Jona-Lasino (ENJL) model 
with chiral U(3) x U(3) symmetry [33-36]. Indeed, all phenomenological low-energy in- 
teractions predicted by Effective Chiral Lagrangians [28,31] for the nonet of low-lying 
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vector mesons (p(770), c<j(780) and so on) can be derived within the ENJL model by 
calculating one-constituent quark loop diagrams at leading order in the long-wavelength 
approximation. As has turned out the long-wavelength approximation works very good 
in spite of the fact that the constituent quark loop diagrams are defined by point-like ver- 
tices of quark-meson interactions and the Green functions of the free constituent quarks 
with constant masses M q ~ 330 MeV, and the masses of vector mesons exceed twice the 
constituent quark mass. A formal justification of the validity of this approximation can 
be given by attracting the Vector Dominance (VD) hypothesis [31,37] due to which the 
effective vertices of low-energy interactions of low-lying vector mesons should be smooth 
functions of squared ^-momenta of interacting mesons varying from on-mass shell to 
zero values. This allows to calculate the vertices of low-energy interactions of vector 
mesons keeping them off-mass shell around zero values of their squared momenta [37], 
and then having had kept the leading terms of the long-wavelength expansion to continue 
the resultant expression on-mass shell. Such a procedure describes perfectly well all phe- 
nomenological vertices of low-energy interactions of low-lying vector mesons predicted by 
Effective Chiral Lagrangians [31,34-37]. 

One cannot say exactly, whether we really have in the RFMD some kind of the VD 
hypothesis, i.e., smooth dependence of effective low-energy interactions of the deuteron 
coupled to other particles on squared 4-momenta of interacting external particles includ- 
ing the deuteron. However, the application of the long-wavelength approximation to the 
computation of one-nucleon loop diagrams leads eventually to effective local Lagrangians 
describing reasonably well a dynamics of strong low-energy nuclear interactions. The 
static parameters of the deuteron and amplitudes of strong low-energy interactions of the 
deuteron coupled to nucleons and other particles can be described in the RFMD in com- 
plete agreement with the philosophy and technique of the derivation of Effective Chiral 
Lagrangians within effective quark models motivated by QCD. 

The agreement between the reaction rates for the neutron-proton radiative capture, 
which is the Ml transition, calculated in the RFMD and the PMA is not surprising. 
Indeed, it is known from particle physics that the radiative decays of pseudoscalar and 
vector mesons like 7r° — > 77, uo — > 7r°7 and so on, caused by the Ml transitions, can 
be computed both in the non-relativistic quark model [38] , which is some kind of the 
PMA, and in the Effective Chiral Lagrangian approach. In the non-relativistic quark 
model the matrix elements of these decays are given in terms of magnetic moments of 
constituent quarks proportional to 1/M q , whereas in the Effective Chiral Lagrangian 
approach they are defined by the axial anomaly and proportional to 1/F n , the inverse 
power of the PCAC constant F n = 92.4 MeV [16,27,39]. Equating the matrix elements 
of these decays calculated in the non-relativistic quark model and in the Effective Chiral 
Lagrangian approach one can express a constituent quark mass in terms of the PCAC 
constant F w [39]. The estimated value of the constituent quark mass M q ~ 400 MeV is 
comparable with the values M q = 330 -j- 380 MeV accepted in the literature [38]. This 
testifies that both the non-relativistic quark model and the Effective Chiral Lagrangian 
approach describe equally well the dynamics of strong low-energy interactions of low-lying 
mesons even if for the decays caused by the Ml transitions. Referring to this example the 
agreement between the reaction rates for the neutron-proton radiative capture calculated 
in the RFMD and in the PMA, respectively, is understandable. The computation of the 
astrophysical factor for the solar proton burning and the disintegration of the deuteron 
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by anti-neutrinos and photons in agreement with the PMA has only confirmed a validity 
of dynamics of strong low-energy nuclear interactions suggested by the RFMD. 

The only problem which is left to discuss concerns the Coulomb repulsion between 
protons in the process of the solar proton burning. As has been shown by Kamionkowski 
and Bahcall [18] the Coulomb repulsion plays an important role for calculation of the 
astrophysical factor S pp (0) = 3.89 x KT 25 MeVb. In our case the Coulomb repulsion is 
taken into account in the form of the Gamow penetration factor C{rj) and, apart from 
weak interactions, only strong low-energy nuclear forces are responsible for the value of the 
astrophysical factor 5 , pp (0) = 4.02 x 1CT 25 MeVb. This discrepancy with the PMA, which 
can be attributed to the peculiarity of the model with a local four-nucleon interaction like 
Eq. (|l.lj ) and a description of strong low-energy nuclear interactions through one-nucleon 
loop exchanges, we are planning to resolve in our further development of the RFMD. 

Finally, we have also shown that in the RFMD with the effective local four-nucleon 
interaction Eq. ( |1 . 1|) one can describe low-energy elastic NN scattering in terms of the 
S-wave scattering length onn and the effective range tnn in spirit of the EFT approach 
[11-13] and in complete agreement with low-energy nuclear phenomenology. 



Appendix A. Computation of the amplitude of the so- 
lar proton burning 

In order to acquaint readers with the machinery of the RFMD we give below the detailed 
derivation of the amplitude of the solar proton burning p + p^D + e + + i/ e . 

The process p + p^D + e + + z/ e runs through the intermediate W-boson exchange, 
i.e., p + p — > D + W + — > D + e + + i/ e . The RFMD defines the transition in terms of 
the following effective interactions [1,2] (see Eq. ( |1.1|) ): 

AipDfX) = -igv[p c {x)l^n(x) - n c (x)^p(x)} D^x), 

1 1-1 1 

{\p(t,^+ -p)Yl 5 p c (t,x--p)} [p c (t,x + -p) 7/ , 7 5 p(t,f - -p)\ 

+\p(t, x + ^p) 7 V (t, x--p)\ [p c (t, x+^p) 7 5 p(t, x - ~ p )]}, 

£ npW (x) = cos$ c [n(x)Y(l ~ 9Al 5 )p(x)} W'{x). (A.l) 

For convenience, in the effective local four-nucleon Lagrangian Eq. (|1 . 1| ) we have intro- 
duced the interaction over a radius-vector p of a relative movement of the protons with 
a 5-f unction 5^(p). 

Then, the transition W + — > e + + v e is defined by the Lagrangian 

A, ee +w(*) = [VV C (^(1 - 7 5 )Va*)] W+(x). (A.2) 
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The electroweak coupling constant <?w is connected with the Fermi weak constant Gp and 
the mass of the W-boson M w through the relation 



9w G 



F 



(A.3) 



In order not to deal with the intermediate coupling constant gw it is convenient to apply to 
the computation of the matrix element of the transition p + p — > D + W + the interaction 

£n P w(*) = [n(x)Y(l - gAl 5 )p(x)] W~(x), (A.4) 

and for the description of the subsequent weak transition W + — > e + + v c to replace the 
operator of the W-boson field by the operator of the leptonic weak current 

W-{x) -> [^(x) 7 ,(l - VHO*)]- (A.5) 
The S matrix describing the transitions like p + p — > D + W + is defined 



J d^x [£ npD (x) + £ npW (x) + ^r pp w + • • •] 



S = Te •/ ' ' en v ' J , (A.6) 

where T is the time-ordering operator and the ellipses denote the contribution of inter- 
actions irrelevant to the computation of the transition p + p — > D + W + . 

For the computation of the transition p + p — > D + W + we have to consider the third 
order term of the S matrix which reads 

S (3) = jj d A Xl d A x 2 d A x 3 T(i£ n M*i) + £n P w(*i) + C^ w ( Xl ) + ...} 

X [C QpB (x 2 ) + £n P w(z 2 ) + £^ PP (X 2 ) + . . .] 
X [C npD {x 3 ) + £npw(x 3 ) + C^T^M + ...]) = 
= -I J d 4 Xl d A X 2 d 4 X 3 T(£ P c r PP (xi)£nM*2)£n P w(x 3 )) + ... (A.7) 

The ellipses denote the terms which do not contribute to the matrix element of the 
transition p + p — > D + W + and the interaction £ np w(x) is given by Eq. (A.4). The S 
matrix element S pp ^ DW+ contributing to the transition p + p — > D + W + we determine 
as follows 

S p 3 pLdw+ = ~i J d A x x d A x 2 d A x 3 T(£ir PP (xi)£n P D(x2)CnM^))- ( A -8) 

For the derivation of the effective Lagrangian £ PP _ >DW + (x) containing only the fields of the 
initial and the final particles we should make all necessary contractions of the operators 
of the proton and the neutron fields. These contractions we denote by the brackets as 

< S pp Ldw+ >= / d A x x d A x 2 d A x 3 < T(£ c P r PP (^i)^np D (a: 2 )/:npw(a:3)) > • (A.9) 

Now the effective Lagrangian £ pp ^ DW + (x) related to the S matrix element < S p _ >DW+ > 
can be defined as 



< S [ 3 p Idw+ >=i J d A xC pp ^ DW +(x) 
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= -ij d 4 x x d 4 x 2 d 4 x 3 < T(Ciy pp ( Xl )C npD (x 2 )C npW (^)) > • (A.10) 

In terms of the operators of the interacting fields the effective Lagrangian £ pp ^dw+( x ) 
reads 

J d 4 x £ PP ^ DW + (x) = - J d 4 x x d 4 x 2 d 4 x 3 < T(£ p ^ pp (xi)£ npD (x 2 )£ npW (x 3 )) > 

= -^G 7rN N x (-igy) x (-g A ) J d 4 x 1 d 4 x 2 d 4 x 3 J d 3 p8 {3 \p) 
x T([p c (t 1 ,f 1 + ±p) la ^ p (t 1 ,x 1 -\p)]Dl(x 2 )W-(x 3 )) 
x < OlTflftfx,:?! + ^p)n 5 P c (h,x 1 -^p)}[pc(x 2 )Yn(x 2 ) - n c (x 2 )^p(x 2 )} 
x [n{x 3 )Yl 5 p{x 3 )})\0 > -^G wNN x (-igy) x (-g A ) J d 4 x x d 4 x 2 d 4 x 3 J d 3 p5 {3 \p) 
x T^tufr + ^ph 5 p(t 1 ,x 1 - l -p)\ D\{x 2 ) W~(x 3 )) 

x < 0|T([p(t 1 ,f 1 + ^p) 1 ^( tl ,x 1 -^p)W{x 2 )^n{x 2 )-n-{x 2 )^p{x 2 )\ 

x[n(x 3 )Yl 5 p(x 3 )})\0>. (A.ll) 

Since p + p — > D + W + is the Gamow-Teller transition, we have taken into account the 
W-boson coupled with the axial nucleon current. 

Due to the relation n c (x 2 )^p(x 2 ) = —p c (x 2 )^n(x 2 ) the r.h.s. of Eq. (A.ll) can be 
simplified as follows 

J d 4 x £ PP ^ DW + (x) = - J d 4 x l d 4 x 2 d 4 x 3 < T(£^" >pp (x 1 )£ n p D (a; 2 )£ n p W (a;3)) > 
= G^nn x {—igy) x g A J d 4 x\d 4 x 2 d 4 x 3 J d 3 p5 ( - 3 \p) 
x Tfl^x, fx + i p ) 7a7 5 p(^i, ^1 - \ p)\ Dl(x 2 ) W-{x 3 )) 
x < OlTflflfx,*! + ip)7 a 7V(*i,xi - ip)][p c (x 2 ) 7 ^(a: 2 )][n(x3)7Vp(^)])|0 > 
+G ffNN x (-ic/v) x #a y rf 4 xirf 4 a;2rf 4 a;3 ^ d 3 p5 (3 \p) 
x T\p(*i,£i + iph^,^ - ip)] £)t(x 2 ) ^-(x 3 )) 

x < 0|T([p(t 1 ,f 1 + ^p) 7 ¥(*i,£i - ip)][p c (x 2 ) 7 ^(a; 2 )][n( a ;3)7 iy 7 5 ^3)])|0 > . (A.12) 
Making the necessary contractions we arrive at the expression 

J d 4 xC pp ^ DW +(x) = - J d 4 xid 4 x 2 d 4 x 3 < T(Cll^ pp (xi)C npD (x 2 )C npAN (x 3 )) > 
= 2 x G wNN x (-igy) x g A J d 4 x 1 d 4 x 2 d 4 x 3 J d 3 p5 {3 \p) 
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x Tfl^f!,*! + 1 -p) lal 5 p(t 1 ,x 1 - l -p)]Dl{x 2 )W-{x,)) 

x (-l)tr{ 7 a 7 5 H)^i ~t2,x 1 -x 2 - l -p)^(-i)S F (x 2 -xz)Yl h 

1 

x (-i)S F (t 3 - t u x 3 - fi - - p)} 
+2 x Gttnn x (-i(? v ) x # A y d 4 Xl d 4 x 2 d 4 x 3 
x TQ^,^ + ±f?) 1 5 p( tl ,x 1 -±p)\ Dl(x 2 ) W~(x 3 )) 
x (-1) tr{ 7 5 H)^(t! - t 2 , x x -x 2 - X - p)^(-i)S F (x 2 - x 3 )Yl 5 

x(-i)SF(t3-ii,i3-^i-^1}, (A.13) 

where the combinatorial factor 2 takes into account the fact that the protons are identical 
particles in the nucleon loop. Let us corroborate the appearance of the factor 2 by a direct 
calculation: 

J d 4 x £ pp ^dw+ (x) = - J d 4 x 1 d 4 x 2 d 4 x 3 < T(Cl^ pp (x 1 )C npD (x 2 )C npW (x 3 )) > 

= G nNN x (-iffy) x 9aJ d 4 x x d 4 x 2 d 4 x 3 J d 3 p5 (3) (p) 
x T([^(t 1 ,f 1 + l -p) lal ?> p (t 1 ,x 1 -±p)]Dl(x 2 )W-(x 3 )) 

x < 0|T([p Q1 (ti,fi + \p) (7 a 7 5 C)ai/3iP/3i(*i>^i -^p)][Pa 2 {x2){C^) a2 fi 2 n f)2 {x 2 )} 
x [n a3 (x 3 )(Yl 5 ) a3 p 3 P(3 3 (x 3 )})\0 > +(7«7 5 ® 7°7 5 -> 7 5 ® 7 5 ) = 
= GttNN x (-i# v ) x # A / d^id 4 ^ 4 ^ / d 3 p5 (3) (p) 

x T([^(t 1 ,f 1 + Ip) 7a7 ^(t 1);?1 - Ip)] /Jj^)^-^)) 

x {(7 a 7 5 C)ai/3i(-0'S'i ;, ( t 2 - *i,x 2 - fi + ^p) a2/3l (C , 7 M ) a2/ 3 2 (-i)S' F (a;2 - XsWs 

x {l v l b ) a ^ 3 {-i)S F (t 3 -t 1 ,x 3 -x l - ]^p)^ ai - 
-(-i)S F (t 2 -h,x 2 -xi - ^p) a2ai {l a l 5 C) ail3l {CY)a 2f 3 2 {-i)S F {x 2 -x 3 )f3 2Ce3 



7 5 <g> 7 a 7 5 -> 7 5 <g> 7 5 ) = 



X (7 !y 7 5 )a 3 /3 3 (-0^( t 3 - fl,&3 - ^1 + T^Wi} + (la 

= G. NN x (-^ v ) x g A J d 4 x ld 4 x 2 d 4 x 3 jd 3 p5^(p) 
x T([p" c (ti,fi + \p) la ^p(t 1 ,x 1 -\p)]Dl(x 2 )W;(x 3 )) 
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x {tr{ 7 a 7 5 CH).Si^2 - h, x 2 - x, + ±p)C>f(-i)S F (x 2 - x 3 )Yl b 

x (-i)S F (t 3 -t!,x 3 -xi- ^p )} 
-(-i)[S F (t2 ~h,S 2 -x, - \p)l a l b C] a2Pl [C^{-i)S F {x 2 -x 3 )Yl 5 
x (-i)S F (t 3 - h,x 3 -x ± + ^p)]a 2 ^} + {lal 5 ® 7°7 5 -> 7 5 ® 7 5 ) = 
= G wNN x (-i^v) x # A f d 4 x 1 d*x 2 d 4 x 3 J d 3 p5 {3) (p) 
x T([^(t 1 ,f 1 + 1^)7^^,^ - ^p)]^,) W-fo)) 
x {(-l)tr{ 7 V[C T ^(t 2 - ti,f 2 -x 1 + ±p)C]YH)S F (x 2 - x 3 )Yl 5 

x (-i)S F (t 3 - h, x 3 - fi - ^p)} 
+(-l)tr{H)[S£(t 2 - *x, £ 2 - fx + ip)7 a 7 5 C] VH)M*2 - x 3 ) 7 "7 5 
x (-i)S>(i 3 - t u x 3 - fx + ^p)}} + ( 7a7 5 <g> 7 a 7 5 -> 7 5 <g) 7 5 ) = 
= G^nn x (-i^v) x # A J d 4 x 1 d 4 x 2 d 4 x 3 J d 3 p5 (3 \p) 
x T\p(*i,£i + ip)7a7 5 p(ti,£i - ^p)]^,) W-(x 3 )) 
x {(-l)tr{ 7 a 7 5 H)^i -£ 2 - I/?)^-^^ -x 3 )Yl b 

x (-i)S F (t 3 - h, x 3 - fi - -p)} 
+ (-l)tr{C T ( 7 a 7 5 ) T CH)^i - t 2 , fi - f 2 + ^p )7 m H)Sf(z 2 - ^hV 

x (-i)S- F (* 3 ~t 1 ,x 3 + x 1 + ^p)}} + ( 7a7 5 ® 7°7 5 - 7 5 ® 7 5 )- (A.14) 

Here we have used the relation C = —C T . Then, by applying the relation C T {j a rf) T C = 
7°7 5 we obtain the following expression 

J d 4 x £ pp ^dw+ (x) = - J d 4 x 1 d 4 x 2 d 4 x 3 < T(Cl^ pp (xi)C npI) (x 2 )C npW (x 3 )) > 
= GWn x (-ig v ) x 9aJ d 4 Xl d 4 x 2 d 4 x 3 J d 3 p5 (3 \p) 
x {Tfl^ix, fx + l - p) 7a 7 5 p(ti,^i - \ P)\ D\{x 2 ) W~{x 3 )) 
x (-l)tr{ 7 a 7 5 H)^i - t 2 , fx - x 2 - l -p)^(-t)S F (x 2 - x 3 )Y^ 
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x (-i)S F (t 3 - h, x 3 - fi - ^p)} 
+T([p-(t 1 , fx + X - p) 7a 7 5 p(ti,^i - \ P)\ Dl(x 2 ) W~(x 3 )) 
x (-l)tr{ 7 a 7 5 H)^i - *2, f i - x 2 + ^p)CY(-i)S F (x 2 - x 3 )Yl 5 

x (-i)S F (t 3 -t 1 ,x 3 + x 1 + ±p)}} + ( 7a7 5 ® 7°7 5 - 7 5 ® 7 5 )- (A.15) 
Using the property of the operators 

[p c (£i, fi + i p ) rp(ti, fi - i p )] = [p c (ti, f i - ^ p) rp(t!, fi + ^ P )] (A.16) 

for T = 7 a 7 5 and 7 5 , we get 

y rf 4 x £ PP ^ DW+ (x) = - J d 4 x 1 d A x 2 d 4 x 3 < T(Cls^ pp (x 1 )jC n pD(x 2 )jC npW (x 3 )) > 

= G-kNN x {~Wv) x 9a J d 4 X\d 4 x 2 d 4 x 3 J d 3 p5 ( - 3 \p) 
x {T([p c (t 1? xi + i p) 7a7 5 p(t!, ^ - i p )] Dl(x 2 ) W'(x 3 )) 
x (-l)tr{ 7 Q 7 5 ^(ti - h,x, -x 2 - l -p)Y{-i)S F {x 2 - x 3 )Yl 5 
x (-i)S F (t 3 - h, x 3 - fi - -p)} 
+T([p c (t 1 ,f 1 - Ip) 7a7 5 p(t 1 ,f 1 + Dj^) ^-(x 3 )) 
x (_l)tr{H) 7 a 7 5 SM*i - h, x 1 -x 2 + ^p)Cr(-i)S F (x 2 - x 3 )Yl* 

x (-i)S F (t 3 -t 1 ,x 3 + x 1 + ^p)}} + ( 7a7 5 ® 7°7 5 - 7 5 ® 7 5 )- (A.17) 

Making a change of variables p — > — p in the last term, we arrive at the expression 
Eq. (A.13). 

Then, S F (x) and S F (x) are the Green functions of the free anti-nucleon and nucleon 
field, respectively: 

r d A k p~ ik ' x 

S F (x) = CS T F {-x)C T = S F (x) = / _- r . (A.18) 

J (Air)' 1 Mn — k 

Passing to the momentum representation of the Green functions we get 

J d 4 x £ PP _ +DW + (x) = 
-_ i04Gl Sv_ /> f d 4 x 2 d 4 k 2 d 4 x 3 d% ik . (a . _ Xl ) jk 3 ■ (x, - x, ) 
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x / dV 3 toT([p _c (ii,£i + ^P)7a7 5 p(*i,£i- lp)]Dl(x 2 )W-(x 3 )) 
x /^i e ^-/ 7 tr{ 7 V- 



d 4 h iff . \ n 5 : 1 _ 5 1 



{ M N -ki + k 2 M N - fcx M N - fci - k 3 J 
igA^NN Sn2 Jd Xl J (2 ^ )4 ^ e 

x I d 3 p^ (3) (p)T([p" c (ti,fi + ^p) 7 5 p(*i,*i -^p)] £>t(a; 2 ) W7(x 3 )) 

x f^e^-PtrL 5 1 — i -1, (A.19) 

J k 2 i { M N -k 1 + k 2 M N - fcj M N - fcj - fc 3 J 

— # — * — * 

where g = fci + (& 3 — k 2 )/2. 

In order to obtain the effective Lagrangian describing the process p + p^D + c + + 
v c we have to replace the operator of the W-boson field by the operator of the leptonic 
weak current Eq. (A. 5): 

yA " NN v^87r 2 i 1 J (2tt) 4 (2vr) 4 

x / d 3 p5 (3) (p)T([p^i,fi + ip)7«7 5 p(*i,*i -\p)\Dl{x 2 ) [^ (x 3 )7,(l -7 5 )^c(^3)]) 
j4 1. . _ r 1 1 1 ~) 



x j^e^-hA^r 

J n z i ! 



_y 7 5 ^_ 

M N - fcx + k 2 M N -hi M N — ki — k 3 



q km—— /d 4 ^ / ^ x 2 d fe 2 d x 3 d A; 3 -jfc 9 • (x, - xi)-ik, ■ (x, - x,) 
+igA^NN V 2 8n2 Jdx 1 J ^ ^ e 

x I d 3 p^ (3) (p)T([p" c (ti,fi + ^)7 5 p(*i^i -\p)]D\{x 2 ) foMxsMl -7 5 )^e(x 3 )]) 

x f^le^-PtrL 5 1 ^7" ^-tV -V (A.20) 

J { M^-k x + k 2 M N - ki M N - fcj - k 3 J 

Now we are able to determine the matrix element of the process p + p — > D + e + + z/ e as 
ci 4 x < D(/c D )e + (/c c +)z/ c (A; 1 ,J|/:p P _ Dc + i , c (x)|p(pi)p(p2) >= 



= (ftr)'fM(fa + *,-»-,») *I(P + P->D + «* + *) _ (A21) 

/2EiF 2£ 2 V 2£ D V 2E e +V2E v y 



where fc^ = k c + + k Uc is the 4-momentum of the leptonic pair, E t (i = 1,2,D, e, u c ) are 
the energies of the protons, the deuteron, positron and neutrino, V is the normalization 
volume. 

Taking the r.h.s. of Eq. (A.20) between the wave functions of the initial |p(pi)p(p 2 ) > 
and the final < D(/c D )e + (/c e +)z/ e (/c I/e )| states we get 

m ar(4)/, \ M (p + P -> D + e+ + U e ) 

(2tt) 5 ( >{k D + k e -p 1 -p 2 ) = 

^2E 1 V2E 2 V 2E D V 2E C+ V 2E Uc V 
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- l9A ^ NN V28^J dXl J (2tt)< (2*)< 6 
x / dV 3) (p) <D(Me + (M^J|T([p-^ 

x [^ e (x 3 )7i/(l - 7 5 )^c(a:3)])b(pi)p(P2) > 
x /^V^tr{ 7 V 1 ^-tV ' ' 



x 



Tr 2 i {' M N -k 1 + k 2 ' M N - ki ' ' M N -ki-k 3 

+ig A G^^ [d* Xl [ d ^ 2d y d ^ d y e~ ik * ■ (*' - *0e-'*3 • (*3 ~ *i) 

x | rf 3 p^ (3) (p) < D(Ato)e + (A^)i/.(^ 

x [Vv e (x 3 )7i/(1 ~ 7 5 )^c(a;3)])b(pi)p(P2) > 

^e^-PtrL 5 1 ^ 1 -V (A.22) 

n 2 i 1 M^-k 1 + k 2 Mn-fa M^-h-fa) 

Between the initial \p(pi)p(p2) > and the final < D(/c D )e + (A; e +)z/ e (A; !/e )| states the matrix 
elements are defined 

< D(A; D )e + (A; e+ )z/ c (A; i , c )|T([^(t 1 ,f 1 + ± p) lal b p{h, x x - X - p)\ (x 2 ) 

X ftMz 3 )7*(l - 1 5 )M X 3)})\P(P1)P(P2) >= [u c {p 2 halMPl)}Hk Uc )lu{l - 1 5 )V{K+)] 

e -i(pi + P2) ■ xi e ik B ■ x 2 e ik e ■ x 3 
x e*(fc D ) V2 ^pp(p) in , 

yj2E 1 V2E 2 V2E B V2E e+ V2E Ue V 

< D(A; D )e + (A; c+ )z/ c (A; 1 ,J|T([^(t 1 ,f 1 + lp)^p{t 1 ,x l - X - p)\ (x 2 ) 

x $v.(x 3 )-y v (l -7 5 )^e(x 3 )])b(pi)p(p 2 ) >= [« c (p 2 hMPi)M^7,(l ~7>(M] 

p -i(Pi + Pt) ■ xt ik D ■ x 2 ik e ■ x 3 
x e;(k B ) V2 ^ p (p)in ; (A.23) 
yj2E 1 V2E 2 V 2E B V2E e+ V 2E Ve V 

where ip pp (p)m is the wave function of the relative movement of the free protons in the 
1 S -state normalized to unit density [40]: 

, . . sin kp 

^ PP (p)in = -j^-, (A.24) 

where k is a 3-momentum of a relative movement of the protons. Since the spatial part 
of the wave function of the protons is symmetric under permutations of the protons, so 
the spinorial part should be antisymmetric. In our approach the spinorial part of the 
wave function of the protons is described by [w c (p2)7<*7 5 «(pi)] and [u c (p 2 )^ 5 u(pi)}, an- 
tisymmetric under permutations of the protons: [u c (p 2 ) r y a 1 5 u(pi)] = — [ uC {pi)lal b u{p 2 )] 
and [w c Go 2 )7Mpi)] = -[u c {pi)l b u(p 2 )\. 



21 



Now let us discuss in details the computation of the matrix elements: 

1 1 

< 0|p c (ti,f! + -p)Fp(t u xi - - p)\p(pi)p(p2) >, 



(A.25) 



where we have denoted T = 7 a 7 5 or 7 s . 

In the quantum field theory approach the wave function \p(pi)p(p2) > should be de- 
scribed in terms of the operators of the creation of the protons a'(pi,<7i) and aT(p 2 ,cr 2 ), 
where pi and Oi {% — 1, 2) are the 3-momenta and the polarizations of the protons. There- 
fore, |p(pi)p(p 2 ) > reads 



\p(Pi)p(P2) >= -^a\p 1 ,a 1 )a ii (p2,a2)\0 > . 



(A.26) 



The wave function Eq. (A.26) is taken in the standard form [41]. It is antisymmetric 
under permutations of the protons due to the anti-commutation relation 

a\pi, <7i) a\p 2 , a 2 ) = ~a\p2, cr 2 ) a f (pi, <7i) 

and normalized to unity. The factor l/y/2 takes into account that the protons are corre- 
lated in the initial state. 

The operators of the proton fields p c (ti, x± + \ p) and p(t±, x\ — \ p ) we represent in 
terms of the plane-wave expansions 



P c (ti,£i + 2 P. 



E 



2E ?1 V 



p(*i,fi - -p) = E 



a(g 2 , a 2 ) «(<&) e"^*! + ' ( f 1 " ^/ 2 ) 



+b\q 2 , a 2 ) e iS **i " ^ 2 ' ( fl " # 2 ) 



(A.27) 



where a(gj, a,) (i = 1, 2) and ai) (i = 1, 2) are the operators of the annihilation and 
the creation of protons and ani-protons, respectively. The computation of the matrix ele- 
ment Eq. (A.25) runs the following way. Holding only the terms containing the operators 
of the annihilation of the protons we get 

< 0|p c (ti,fi + ^p )Tp(t 1 ,x l - ^ p)\p(pi)p(p2) >= 



<?l,"l 9*2, "2 



1 



i p ~Kqi + 92) • zi + - <f 2 ) • p/2 



v <?2 



1 



x [w c (gi) r«(g 2 )] < 0|a(gi,ai)a(g 2 ,a 2 )a t (pi, 0-1)^(^2,0-2)10 > 



(A.28) 
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e 



The vacuum expectation value < 0|a(g 1; a±) a(gi, ai) a\f>i, <Ji) o){p 2 , cr 2 ) |0 > reads: 

< 0\a(qi,a 1 )a(q 1 ,a 1 )a^(p 1 ,a 1 )a j (p 2 ,a 2 )\0 >= 

$qipi $ct\a\ $q2P2 ^a20"2 ^<?2Pi ^a2ci ^qiP2 ^aio"2) (A. 29) 

where we have used the anti-commutation relations 

a(q, a) a^(p, a) + a\p, c)a(?, a) = S^p5 aa (A. 30) 

and the properties of the operators of the creation and the annihilation: < 0\a^(p, a) = 
and a(q, a)\0 >= 0. 

Substituting Eq. (A. 29) in Eq. (A. 28) and summing up the momenta and the spinorial 
indices we arrive at the expression 

1 i e ~i(pi + P2) ■ xi 
< 0|p c (ti,fi + -p)r>(£i,£i - -p)\p{pi)p{p2) >= * 

2 2 ^2E 1 V2E 2 V 

x -L (^uc{ pi )Tu{p 2 )]e i ^-^)-pl' 2 - [u^Tuip^e- 1 ^ ~ ' p/ 2 ^ = 

. v^2 [u c (p 2 ) T u( Pl )] - Bigg{e % ^ - P2) • + e" l (Pi - ^) " P/2 , 

^2E 1 V2E 2 V 2 ) 

(A.31) 

where the relation [w c (pi) T -u(p 2 )] = — [-u c (p 2 ) T -u(pi)] has been used. The sum of the 
exponentials 

I ^(Pi - P2) ■ p/2 + e -i(pi - p 2 ) ■ p/2^ 

describes the spatial part of the wave function of the relative movement of the free protons. 
This wave function is symmetric under permutations of the protons and normalized to 
unit density [40]. Since the protons should be in the 1 So-state, expanding exponentials 
into spherical harmonics and keeping only the S-wave contribution we obtain [40]: 

Ujk-? +e -ik.fi\ = ™kp + ^ (A . 33) 
2 V J kp 

where k = (pi — P2)/ 2 is the relative momentum of the protons. This completes the 
explanation of the derivation of the matrix elements in Eq. (A. 23). 

Substituting the matrix elements Eq. (A. 23) in the r.h.s. of Eq. (A. 22) we obtain the 
matrix element of the solar proton burning in the following form 

(2ir)^\k D + k e - Pl - p 2 ) iM(p + p -> D + e + + v e ) = 
G 

= -V2C(rj) ffA^NN - J 2 [M c (p 2 )7 a 7 5 «(Pi)] [u(k Vc )j v (l - -f 5 )v(k c +)} e* (k D ) 



x f d 4 xi ( - X2d k<2 - Xsd ks e i( ^ 2 + k 3~Pi- P2) ■ xi J(k D - k 2 ) ■ x 2 i(k e - k 3 ) ■ x 3 
J J (2nY (2tt) 4 
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^sinfcp [dkijq.pj^, 1 ^ ^Yl 5 - 



x / £p8^(pf-^ ( Ip-eW-Ptr 7 V 
J fcp J it i { M N -k 1 + k 2 M N - ki M N - fci - fc 3 , 

- v ^C'(7 ? )^ A aNN^^[w c (j9 2 )7 5 w(pi)][«(A;,j7,(l-7 5 )^(A: c+ )]e;(A; D ^ 

x /" ^ f dAx 2d 4 k2 d 4 x 3 d 4 k 3 ^ + kz _ pi _ p2 ) . Xl ^ _ k ^ . X2 _ k ^ . Xz 
J J (2tt) 4 (2tt) 4 



sinkp fd% jff.jj^ 1 ^ ^ 1 ^ ^ 5 1_ 



fcp •/ 7r 2 i [ Mn — k\ + k 2 M N - fci M N - fci - fc 3 J 

(A.35) 

where we have appended the Gamow penetration factor C(^) taking into account the 
Coulomb repulsion between the protons [2]. 

Integrating over x±, x 2 , x 3 , k 2 and fc 3 we obtain in the r.h.s. of Eq. (A.35) the 8- 
function describing the 4-momentum conservation. Then, the matrix element of the p + 
p — > D + c + + u c process becomes equal 

iM(p + p -> D + e + + u e ) = -V2C(rj) gAG wNN -^=-^ 

x [u c (p 2 ) lal b u( Pl )}[u(k Vc ) lu (l - 7 >(fc e+ )] e ;(fc D ) 

fd 4 k, f a 5 1 M 1 „ 5 1 1 

x /^-tW 7 V s ^7^ ^7 7 * ^} 

tHz 1 Mn — k\ + k 2 M N - fci M N - fci - fc 3 J 

-V^C'(r ? )^NN^^[« c (p 2 )7 5 «(pi)][«(fc, c )7,(l-7 5 )^(fce + )]e;(fc D ) 

x / ^tr{ 7 5 ! ^7"^7"7 5 i A, (A.36) 

•/ tt 2 z 1 M N - fci + fc D M N - fci M N - fci - fcj 

where we have integrated over a relative radius-vector p too. It is convenient to represent 
the matrix element Eq. (A.35) in terms of the structure functions J afiU (ki),ke;Q) and 
J^(k D ,k e ;Q): 

iM{p + p -> D + e + + u e ) = 



= - C( V ) G Y g A G^f^[u c (p 2 h a l 5 u(Pi)}WK c hv(l ~ l>(k c+ )} e ;(fc D ) J a ^(k D , k e ; Q) 

C( V ) G v g A G wm ^[u c (p 2 hMPi)]WK hv(l - 7 5 )^(fce + )] e;(fc D )^(fc D , k t ; Q), 

(A.37) 



where the structure functions J a ^(kj), ke; Q) and J^fco, kf, Q) are defined as [2] 

J a ^(k D ,k e] Q) = 



fd 4 k f 



1 1 1 



fd 4 k r 



M N -fc-0 + fc D M N -fc-Q M^-k-Q -k e 
J^(k D ,k e ;Q) = 



7 5 ^— : ^7" 1 ^— ^ ^1. (A.38) 

M N -k-Q + k D M N — k — Q M N -k-Q-k e j 
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We have introduced a 4-vector Q = a k^ + b ke caused by an arbitrary shift of a virtual 
momentum with arbitrary parameters a and b. 

Thus, the problem of the computation of the matrix element of the p + p — > D + 
c + + u e process reduces to the problem of the computation of the structure functions 
Eq. (A. 38). Since the energy of the leptonic pair is small compared with the nucleon 
mass, we can set in the integrand kg = [2]. This gives 

J a ^(k D ,k e ;Q) = 

^tr/ 7 V - 7^ - 7 V - 1 

vr 2 i ( M N -k-Q + k D M N -k-Q M N — k — Q ) 

J^(k D ,k e ;Q) = 

\ 4, (A.39) 

M N -k-Q + k D M N — k — Q M^-k-Q] 

For the calculation of the momentum integrals we would follow the philosophy of the 
derivation of Effective Chiral Lagrangians within effective quark models motivated by 
QCD [33-36], in particularly, Chiral perturbation theory at the quark level (CHPT) g [35] 
formulated on the basis of the ENJL model induced by the effective low-energy QCD 
with linearly rising confinement potential [42]. In (CHPT) g all low-energy vertices of 
meson interactions are determined by constituent quark loop diagrams with point-like 
quark-meson vertices and the Green functions of the free constituent quarks with con- 
stant masses M q = 330 MeV [35]. To the computation of the momentum integrals one 
applies a generalized hypothesis of Vector Dominance [31,37] postulating a smooth depen- 
dence of low-energy vertices of hadron interactions on squared 4-momenta of interacting 
mesons. Due to this hypothesis one can hold all external particles off-mass shell at 
squared 4-momenta p 2 much less than M 2 , i.e., M 2 ^> p 2 . Then, after the computation 
of the momentum integrals at leading order in long-wavelength expansion, i.e., in pow- 
ers of external momenta, the resultant expression should be continued on-mass shell of 
interacting particles. Within the framework of this procedure one can restore completely 
all variety of phenomenological vertices of low-energy meson interactions predicted by 
Effective Chiral Lagrangians [28,31,34-36]. It is important to emphasize that this pro- 
cedure works good not only for light mesons like 7r-meson, which mass is less than the 
mass of constituent quarks, but for vector mesons like p(770), u(780) and so on, which 
masses are twice larger than the constituent quark mass. Since the former resembles the 
RFMD, where the mass of the deuteron amounts to twice the mass of virtual nucleons, 
we expect that the long-wavelength approximation should work in the RFMD as well as 
in the effective quark models with chiral U(3) x U (3) symmetry applied to the derivation 
of Effective Chiral Lagrangians. 

Thus, for the computation of the momentum integrals we assume that the deuteron is 
off-mass shell and Mn ^> \J~k% y Then, we expand the integrand of the structure functions 
Eq. (A.39) in powers of k^> keeping only the leading contributions. The result of the 
computation we continue on-mass shell of the deuteron k\^ — > [2]. 

Keeping the leading terms of the expansion in powers of k® we get [2]: 

J a ^(k D , kg; Q) = 3 W - k^) + 1(1 + 2a) W + k^ a ), 
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J^(k D ,k £ ;Q) =jT4M N J 2 (M N ), (A.40) 

where due to the relation &d • e D (&o) = the terms proportional to have been dropped 
out. Then, ^(Mn) is a logarithmically divergent integral defined in the RFMD in terms 
of the cut-off A D = 68.452 MeV such as A D < M N [1,2]: 



f d 4 k 1 F d\k\k* 2 I A D \ 

^ = I ^(M^W = 2 / MTI^ - 3 U) « (A ' 41 > 

The cut-off Ad restricts 3-momenta of the virtual nucleon fluctuations forming the phys- 
ical deuteron [1,2]. Due to the uncertainty relation Ar Ad > 1/2 the spatial region of 
virtual nucleon fluctuations forming the physical deuteron is defined by Ar > 1.44 fm. 
This agrees with the range of nuclear forces (NF) caused by the one-pion exchange with 
the mass M n = 135 MeV: r NF = 1/M n = 1.46 fm [14]. 

After the continuation of the results of the calculation of the structure functions 
on-mass shell of the deuteron the contribution of J^ v (k^ik^Q) can be neglected rel- 
ative to the contribution of J afiU (k D ,kf,Q). The contribution of the structure function 
J afil/ (k D , kf, Q) does not depend on the mass of virtual nucleons and according to Ref. [17] 
can be valued as the anomaly of the AAV one-nucleon triangle diagram. The ambiguity 
of the calculation of ^T afiu (k-D, ke; Q) caused by the dependence on an arbitrary shift of 
a virtual momentum can be fixed by requirement of gauge invariance of the amplitude 
of the process p + p^D + e + + z/ c under gauge transformations of the deuteron field 
e* M (/c D ) — > e* fl (k D ) + A/c£), where A is an arbitrary parameter. This gives a = —1/2 and 
the structure function in the form [2]: 

J "" (km, h- Q) = 3 W - k^gn- (A.42) 

The attraction of requirement of gauge invariance in order to remove ambiguities of the 
structure function J a ^ v {k-£,, ke; Q) and to fix the contribution of the anomaly is in complete 
agreement with the derivation of the Adler-Bell-Jackiw axial anomaly performed in terms 
of one-fermion loop diagrams [16]. 

Since we strive to draw a similarity between the RFMD and effective quark models 
motivated by QCD applied to the derivation of Effective Chiral Lagrangians, requirement 
of gauge invariance under gauge transformations of the deuteron field 

D^x)^D^x) + dJ(x), (A.43) 

where f(x) is a gauge function, can be justified by referring to a dynamics of vector meson 
fields in these effective quark models [34-36]. The effective Lagrangian of the physical 
deuteron field D^x), which we apply to the calculation of one-nucleon loop diagrams 
describing effective low-energy interactions of the deuteron coupled to nucleons and other 
particles, reads [1,2]: 

C(x) = ~ Dl(x)D^(x) + M^Dl(x)D^x) - ig Y \p(x)Yn c (x) - n(x)^p c (x)} D^x) 

—igy[p c (x)^n(x) —n c (x)^p(x)} DUx) +p(x)(i^d IJi — M N )p(x) + n(x)(i^d^ — M N )n(x), 

(A.44) 
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where D^ v {x) = d^D u (x) — d u D fl (x). Since the deuteron field D^{x) couples to the con- 
served nucleon current J M (x) = — igy [fin(x)^rf(x) — n(x) r y tJ 'p c (x)}, i.e., d^J^{x) = 0, 
invariance of the effective Lagrangian Eq. (A. 44) under gauge transformations of the 
deuteron field Eq.(A.43) is violated only by the mass term. The same problem encounters 
itself for description of dynamics of vector meson fields (p(770), c<j(780) and so on) in 
effective quark models with chiral U(3) x U(3) symmetry [34-36]. Since these mesons 
are massive, the kinetic Lagrangians of vector meson fields are not invariant under gauge 
transformations of these fields. Nevertheless, for the derivation of effective low-energy 
interactions of vector mesons coupled to other particles requirement of gauge invariance 
turns out to be very important. For example, by virtue of requirement of gauge invariance 
under gauge transformations of p(770) and c<j(780) meson fields one can fix unambigu- 
ously the coupling constant of the upir interaction defined by the Adler-Bell-Jackiw axial 
anomaly [16] that plays an important role for the correct description of the u — > 7r + 7r~7r° 
decay and many other low-energy processes. 

Substituting the structure function Eq. (A. 42) in Eq. (A. 37) we obtain the matrix 
element of the solar proton burning in the following relativistically invariant form 

iM(p + p^D + e+ + u e ) = -C{r]) G v g A G nNN ^ 

x (AftT - feSiT)[« c (p 2 )7«7 5 «(Pi)][fi(^)7 v (l " 7>(M1 e;(k D ). (A.45) 
In the low-energy limit due to the low-energy reduction 

[u c (p 2 hal 5 u(Pi)} -> -9ao[u c (p 2 h 5 u(Pi)} (A.46) 
the matrix element of the solar proton burning can be brought up to the form 

iM(p + p -> D + e + + v e ) = C(rj) G v g A M N G nNN ^ 

x [u c (p 2 ) i b u(Vi)mh c )^(l - 7 XM] e^ko), (A.47) 

where we have set = Md ~ 2M^ valid on-mass shell of the deuteron. This completes 
the calculation of the matrix element of the solar proton burning. 

Omitting the Gamow penetration factor C(r)) the residual part of the matrix element 
of the solar proton burning Eq. (A.45) can be defined by the effective Lagrangian 

C pp ^ Dc+Vc (x) = g A G^^^ Dl(x) fr c (x)Yl 5 p(x)} [ieW/(l " 7 5 )^e(x)]. (A.48) 

This Lagrangian is local in accordance with the condition of microscopic causality [32] . 

The calculation of the matrix element of the disintegration of the deuteron by anti- 
neutrinos z/ c + D^e + + n + n can be carried out by an analogous way and defined by 
the same structure functions. The effective Lagrangian describing the matrix element of 
the transition z/ e + D^e + +n + n can be written as follows 

£p e D^ e+ nn(z) = SaG»nn^ D^x) [^/^(a;)] [^(^"(l - 7 5 )^c(^)]. (A.49) 

V2 oir z 

The effective Lagrangians Eq. (A.48) and Eq. (A.49) testify distinctly that the processes 
of the solar proton burning and the disintegration of the deuteron by anti-neutrinos are 
governed by the same dynamics of strong low-energy nuclear interactions in agreement 
with charge independence of the weak interaction strength. 
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Appendix B. Low— energy elastic NN scattering in the 
RFMD 



In this Appendix we show how in the RFMD one can describe a phenomenological am- 
plitude of elastic NN scattering by using the effective four-nucleon interaction Eq. 
For simplicity we suggest to consider the elastic low-energy np scattering. In the RFMD 
the amplitude of the elastic low-energy np scattering can be written as follows 

M(n V - np)(A:) = -A(k) — [«(p' 2 )7 5 w c (K)] [u c (pi)l 5 u(P2)}, (B.l) 

since at low energies the interaction ® 7 M 7 5 reduces to 7 5 <8> 7 5 , then pi and p\ 
(i=l,2) are 4-momenta of the proton and the neutron in the initial and final states and 
A; is a relative 3-momentum of the np system. The phenomenological amplitude of the 
low-energy elastic np scattering A(k) ph reads [12,14] 

A(k) ph = — p , (B.2) 

7 2 • 7 

r ~ Triple 
CL np A 

where a np and r np = (2.75 ± 0.05) fm are the S-wave scattering length and the effective 
range of the np-scattering in the 1 S -state [10]. At k — ► we get *4(0) p h = — a np which 
gives the cross section equal cr (np — > np) = 47ra„ p . 
In the RFMD due to the low-energy reduction 

[n(x)^ 5 p c (x)\ [p c (x)7 Q 7 5 n(x)] -> -[n(x)^ 5 p c (x)\ [p c (x) ^n(x)} (B.3) 

the np scattering runs through the one-nucleon loop exchange. Using the effective inter- 
action Eq. fll.lD we can write down the effective Lagrangian for the low-energy elastic np 
scattering: 

J dxL eS (x) scattcrmg - — J d x J (2?r)4 e 

x {[n(x)7 Q7 y(x)][p c (x 1 )7^7 5 n(x 1 )] l 7 Q/3 (A; 1 ) + [n(x)j 5 p c (x)]\p e (x 1 )'y s n(xi)] J(h) 
+ Kx)7 Q 7¥(^)M^ (B.4) 

where J afi (k\), J{k\)-, J a {k\) and J a {k\) are the structure functions defined by the 
momentum integrals: 

J aP {h) = J ^tr{ TT ^- 7 a 7 5 ^ \ ~ 7V^ 



H [M N -q M N -q-kx 
— • 1 M N - q Mn — q — k\ 



28 



■kH [M N -q M^-q-ki ) 

The amplitude of the low-energy elastic np scattering defined by the effective Lagrangian 
Eq. (B.4) reads 

Al(np^np) = -^f 

x {[«0/ 2 )7a7 5 « c (pi)] [u c (pih^ 5 u(P2)] J a \P) + K^hV^i)] [u^ Pl h 5 u(P2)] J{P) 

(B.6) 



where P = p\ + p 2 = p[ + p' 2 and in the center of mass frame P M = (2yk 2 + M^, ). 

Due to the low-energy reduction Eq. (B.3) the amplitude Eq. (B.6) reduces to the 
form 

A<(np - np) = [u(p' 2 h 5 u%p[)} [u^p^u^)} 

x [-J°°(P) + J{P) + J°(P) - J°(P)}, (B.7) 

where the structure functions J 00 (P), J(P), J°(P) and J°(P) are given by Eq. (B.5) 
with the change k± — > P. 

The integrals over q are both quadratically and logarithmically divergent. In the 
RFMD they are regularized by a cut-off A D <C M N [2]. Then, calculating the integrals 
over q we have to take into account that quadratically divergent integrals regularized by 
a cut-off are defined ambiguously with respect to the shift of virtual momenta. Indeed, 
it is well known [17,2] that 



/ 



^trj !— T V !-^7V' 

nh (M„-q-Q M N - q - Q - P 



^tr(— — :7 V - \ 1 I>A + 2[Q-(Q + Pf + QI>(Q + Pr-Q.(Q+P)g^, 
7r z ? [ M N — q M N — q — P J 

dV f 1 5 1 



, . tW ^7 5 s ^7 5 

vr 2 2 \ M N — q — Q M N -q-Q-P 



^tr(,^ 7 ^^--^-2Q.(Q, P). (B.> 



y 5 - ~ 5 

"-/' [ M N - g M N -q-P 

where a 4-momentum Q defines an arbitrary shift of a virtual momentum q — > q + Q. The 
most general form of Q reads : Q = £ P + iV, where £ is an arbitrary parameter and N 
is an arbitrary 4-momentum orthogonal to P, i.e., P ■ N = 0. As P is a time-like vector, 
P 2 > 0, so A" is a space-like, A^ 2 < 0. In the center of mass frame we have = (0, N ). 
The result of the calculation of the structure function — J 00 (P) + J{P) can be given in 
the following general form: 

[- J 00 (P) + .7(F)] = ^(d + C 2 A; 2 ), (B.9) 

where Cj («=1,2) are arbitrary constants containing all uncertainties induced by shifts of 
virtual momenta. One should take into account that virtual momenta in the momentum 
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integrals defining the structure functions J~ al3 (P) and J{P) can be shifted independently. 
The structure functions J7"°(P) and S°{P) are logarithmically divergent and, therefore, 
do not depend on the shift of virtual momentum. The contribution of these structure 
functions can be absorbed by the constants C\ and C 2 . 

9hm [_ JOO(P) + J( P ) + J\ P ) - J\ P )\ = ^-{C 1 + C 2 k 2 ), (B.10) 

Inserting Eq. (B.10) in the r.h.s. of Eq. (B.7) we obtain the amplitude of the low-energy 
elastic np scattering in terms of the constants C\ and C 2 : 

4<7T 

M(np - np) = - vpCCi + C 2 k 2 ) [u(p' 2 )t V ! (p£)] [u^p^u^)}. (B.ll) 

M N 

Neglecting the terms of order 0(k 4 ) as it is accepted for the description of low-energy 
elastic np scattering [12], we arrive at the expression 

4n 1 

M(np np) = -— r [u(P2h 5 u c (p'i)} [u c ( Pl )^u(p 2 )}. (B.12) 

Matching Eq. (B.12) with Eq. (B.l) we obtain ^4(A;)rfmd in the form: 



1 



-4(A;)rfmd = -j r< • (B.13) 

' ° 2 k 2 



Following [12] we can set 

C\ = -a np , C 2 = -^r np al p . (B.14) 
This yields the amplitude of the low-energy elastic np scattering in the form 

*4(*0rfmd = — 1 — \ • (B.15) 

"^ + 2 rnpfc2 

Up to the imaginary part — % k which can be appended to the r.h.s. of Eq. (B.15) due to uni- 
tarity the amplitude -4(/c)rfmd coincides with the phenomenological amplitude Eq. (B.2). 
Thus, we have shown that in the RFMD with the local four-nucleon interaction Eq. (|1 . 1| ) 
one can describe, in spirit of the EFT approach [11-13] and in agreement with low-energy 
nuclear phenomenology, the amplitude of the low-energy elastic np scattering in terms 
of the S-wave scattering length a np and the effective range r np . This confutes the state- 
ment by Bahcall and Kamionkowski [43] that the effective local four-nucleon interaction 
Eq. (|1.1|) leads in the RFMD to the zero effective range for elastic NN scattering, i.e., 
t nn = 0. This statement [43] is also confuted by our results on the astrophysical factor 
for the solar proton burning, the reaction rate for the neutron-proton radiative capture 
and the cross sections for the low-energy disintegration of the deuteron by photons and 
anti-neutrinos obtained in good agreement with the PMA. 
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